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Abstract 

Neutrino mean free paths in magnetized neutron matter are calculated using the Hartree-Fock approxima- 
tion with effective Skyrme and Gogny forces in the framework of the Landau Fermi Liquid Theory. It is shown 
that describing nuclear interaction with Skyrme forces and for magnetic held strengths logio B{G) ^17, the 
neutrino mean free paths stay almost unchanged at intermediate densities but they largely increase at high 
densities when they are compared to the field-free case results. However the description with Gogny forces 
differs from the previous and mean free paths stay almonst unchanged or decrease at densities [1 — 2]po. 
This different behaviour can be explained due to the combination of common mild variation of the Landau 
parameters with both types of forces and the values of the nucleon effective mass and induced magnetization 
of matter under presence of a strong magnetic field as described with the two parametrizations of the nuclear 
interaction. 

PACS numbers: 21.30.Fe,21.65.Cd,26.60.-c. 



I. INTRODUCTION 



The study of the behaviour of hadronic matter in the density-temperature, (p,T), diagram 
allows to have a deeper understanding of matter under extreme conditions. In this context, the 
high density, low temperature limit can be addressed for a fermion system using the Landau Fermi 
Liquid Theory (FLT) From a theoretical point of view the properties of this type of normal 
quantum systems can be studied calculating the interaction matrix element of quasiparticle (qp) 
excitations close to the Fermi surface. The inclusion of an additional component in the problem, 
a magnetic field, B, allows further testing the properties of magnetized Fermi Liquids. The role 
of magnetic fields in bulk properties and equation of state has been partially analyzed in the past 
for nuclear matter and quark matter [J] [5J . Due to the tiny value of the neutron magnetic 

moment (M n = -1.9130427(5)^ (a*jv = 3.1524512326(45) xl(T 18 MeV G" 1 ) [6] and in order to 
provide a sizable magnetization, huge magnetic fields are needed. 

The only scenarios where we have indication of such intense fields are, first, from estimates of 
the background magnetic fields created in heavy-ion collisions like those at RHIC [7] and, second, 
in a subgroup of pulsars called magnetars. For these astrophysical objects surface magnetic field 



strengths are of the order B ~ 10 15 G Q, Is]. Recent numerical simulations [10] of formation of 
proto-neutron stars show that the field configuration plays a significant role in the dynamics of the 
core if the initial magnetic field is large enough. In particular, in the rapid cooling of the newly 
formed neutron-rich object neutrino transport is an important ingredient JjJ]. However, some 



of these simulations lack from accurate and consistent neutrino transport, missing the impact of 
magnetic fields in the microphysics input that affects the dynamics of the collapsing dense objects. 

In most of the existing calculations of nuclear matter (either symmetric, pure neutron or beta- 
equilibrated) the effect due to the presence of strong magnetic fields and the consistently induced 
spin polarization are discarded in a first approximation. Either relativistic [12I \\\ or effective 
approaches have been used to obtain some insight into the equation of state (EOS) or some 
structure properties [2] in presence of magnetic fields. These include a possible transition to a fer- 



romagnetic state, although simulations using realistic potentials seem to 



a non-vanishing magnetization in a low temperature nuclear plasma 
some degenerated observables as obtained in the context of the FLT 



jO prevent it [15(. In general, 
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II. FORMALISM 



In this work we are interested in the response of a spin-polarized pure neutron plasma to a 
weak neutrino probe. It can be seen 3] that for the density range p < 4po, where the quark 
deconfinement is not expected to take place, and for magnetic field strengths of maximum strength 
B ~ 10 18 G, allowed in principle by the scalar virial theorem, the neutral system is mostly neutrons. 
The maximum magnetic field strength we will consider is B* ~ 2 x 10 4 (as measured in units of the 
electron critical field B* = BjB c e with B c e = 4.4 x 10 13 G) and the neutron fraction is Y n > 0.98 



121 ] . So the neutral plasma is mostly neutrons but leptons and additional baryons are also present 
in a tiny fraction that should be considered for full application in an astrophysical scenario where 
/3-equilibrium holds. 

We are interested in exploring the effect of a strong magnetic field and the spin polarization of 
a pure neutron plasma through the structure functions, which provide information on density and 
spin density in-medium correlations. The homogeneous system under study is under the presence 
of an internal magnetic field, B = populated by species with paricle density p a , where a = ±1 
is the spin z-projection. A = p+ ~ p ~ is the spin excess and p = p + + p_ is the total particle density. 
For given thermodynamical conditions A is obtained by minimizing the Helmholtz free energy 
per particle, f(p,T,B,A) = e — p n ApB, where e is the energy per particle. Note that parallel 
(antiparallel) aligned magnetic moments (spins) are energetically favoured. We have considered an 
effective approach to describe the nuclear interaction using zero-range Skyrme forces 
of the most widely used parametrizations given by the Lyon group SLy4 and SLy7 
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2J and 



finite range Gogny with DIP [22j and D1S 23J forces. All of them provide good values for binding 
of nuclei and also for neutron matter EOS. 

In the context of the FLT the properties of non-magnetized systems at low temperature have 
been evaluated 



241 ] by calculating the qp matrix element around the Fermi surface where the 
only dependence is on fermionic densities and the qp scattering angle, 6, involved. In the usual 
formalism, for the non-magnetized case the qp matrix element is written as a multipolar expansion 
in Legendre polinomials, 



V ph = J2lfl + 9i<ri-<r2]Pi(cos6), (1) 

1=0 

fi and gi are the so-called Landau parameters of multipolarity /. In the more general case where any 
two possible spin orientations (a, a') are taken into account, the polarized qp matrix elements 25] 
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[17( are a crucial ingredient to compute the response functions to a a weakly interacting neutrino 
probe that excites a collective mode (uj, q) under the presence of a magnetic field B. The Lindhard 
function in the system, x^ a ' a ^{ UJ }l)j satisfies the Bethe-Salpeter equation and can be written in 
the dipolar (I < 1) case in the random phase approximation (RPA) as a coupled system, 

a"=+- 

+7 w E /Mr^".^), (2) 

o"=+,~ 

.%^U>) V f(™") v (v"v') 



+7r e /r CTj r (CT CT) , (3) 

(T"=+,- 

with the auxiliar definitions, r( CT > CT ') = / j0pcos(9) G ( - a ' a '^ and 7 1 ct) = f Acos n (0) G^ . Notice 
that the qp pro pag ators G^ have been given in [2g] and the expressions for the coefficients 
can be written [25] in the Landau limit as 7^ = v^Xo an d 72°^ = ^^Xo ~~ ^eff 2 where 
^ = T^T%\ ^ ne IP effective mass in a magnetized system depends on the polarized dipolar 
coefficients [l|], 

K/m = 1 + ~N 0a [fi a ' a) + &^)f[^ ] ] (4) 

where NQ a = m ^^ , ' T is the quasiparticle level density at each polarized Fermi surface with momen- 
tum kp,(T- 

The generalized parameters 3 are obtained by derivating the Helmhotz free energy with 



respect to the polarized density component, /k,o-.kV = g nk 9 dn y ., fl?]- setting momenta on the 
polarized Fermi surfaces and expanding the resulting expression as a series in Legendre polinomials 
of multipolarity I. These generalized parameters fullfill the following relations recovering the usual 
ones in FLT in the limit A — > 

ft = ^-^f , (5) 

m = * f • (6) 

With the generalized paramters and using the expressions in Eq. §2$ the corresponding Lindhard 
function for the isovector (S = 0) response of the plasma can be written as, 

x (5=o) = + x (-) + x (+-) + x (-+) f (7) 



and for the vector-axial (S = 1) response as, 

x (5=l) =x ( ++)+x (-)_ x ( + -)_ x (- + ) (g) 

rn 

Then the previous expression of the Lindhard function in RPA approximation [271 ] include in- 
medium correlations at zero temperature. From them, one can obtain the structure functions 
given by, 

S s =°> 1 (^q) = —Irn X S=0 ' 1 (u;, q ). (9) 

7T 

The structure function allows to calculate the non-relativistic differential cross section of neutrinos 
scattering off matter via neutral currents from 2S ] 



= S^ ,2[ ^ (1 + cos9 ) S °^>9) + CK 3 " cosO^icq)} (10) 

where E(E') is the incoming (outgoing) neutrino energy and k (k f ) is the neutrino incoming (outgo- 
ing) three-momentum. The transferred energy is uj = E — E' and the transferred three-momentum 
is q = k — k'. The neutral current vector and axial vector charges are Cy = 1/2 and Ca = — 5a/2 
where g a = 1.260 [6]. Gf/(^c) 3 = 1.166 39(1) x 10~ 5 GeV~ 2 is the Fermi coupling constant. Once 
the response has been evaluated it is straightforward to evaluate the neutrino mean free paths in 
the medium , A -1 = / y J^ dVLduj. 



III. RESULTS 



In this section we include the effect of in-medium correlations in the neutron magnetized system 
as obtained in the Hartree-Fock approximation in the presence of a strong magnetic field. In Fig. [1] 
the ratio of effective neutron mass as compared to the free value at saturation density, po, is shown 
as a function of the logarithm of the magnetic field strength for the Skyrme SLy7 (a) and Gogny 
DIP (b) parametrizations. For each model upper (lower) curves correspond to spin up (down) 
polarized particles. With Skyrme description the intense field affects more both effective nucleon 
mass absolute and relative (up-down polarized components) values than with Gogny. The impact 
of density effects on the mean free path can be seen in Fig. [2j We plot the ratio of change of 
neutrino mean free paths in the RPA dipolar approximation for a fixed value of magnetic field 
strength, B = 5 x 10 17 G, with respect to the field-free case, Rb = A - B=5x1 ^ 7 ^ ^ ag a f unc tion 
of the density (in units of nuclear saturation density, po). We consider Skyrme SLy7 (solid line), 
SLy4 (long dashed line) Gogny DIP (short dashed line) and D1S (dotted line) parametrizations 
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and set as a typical value of neutrino incoming energy E v = 15 MeV. While Gogny forces show 
almost unchanged or very mild reduction of mean free paths at densities in the range [1 — 2]po, the 
Skyrme forces show a high density dramatic increase with respect to the field-free case. Note that 
all standard Skyrme forces predict the onset of a ferromagnetic transition in the range [1 — <i]po, 
and in our selection of interactions for the study cases it is near 3.3/?o- However this feature is 
not present in the Gogny forces that prevent ferromagnetic transitions. At densities [1 — 2]po, 
effects due to the energetic contribution of the magnetic perturbation, introduced by the neutron 
magnetic moment, are small at the selected field (B = 5 x 10 17 G) with respect to changes in other 



single particle 
(see Fig. 7 in 



properties like effective masses. For even lower densities (i.e. 0.5po) it can be seen 



lq) that the /i n ApB term produces a relevant contribution to the magnetization. 



For application in astrophysical scenarios and at low densities one should consider that the effect of 
a non-zero proton fraction determines the appearance of pasta phases (29 j where electromagnetic 
and nuclear interactions are frustrated and clustering of matter arises. As density grows, at fixed 
values of B, the spin polarization decreases forming a plateau at intermediate densities before the 
possible appearance of a phase transition in the system. 

In Fig. [3] we plot the ratio of change of neutrino mean free paths in the RPA dipolar approxima- 
tion computed for a generic value of magnetic field strength with respect to the field-free case as a 
function of the logarithm (base 10) of the magnetic field strength, R p = ^ B \g B Q =0 ■ We set a value of 
density p = 3po and use SLy7 (solid line), SLy4 (long dashed line), DIP (short dashed line) and D1S 
(dotted line) parametrizations. For fields below B ~ 10 17 G there is almost no change in the ratio 
but for larger strengths there is a decrease (increase) as computed with Gogny (Skyrme) forces. 
For this high density case the change can be ~ 10% as computed with the SLy7 parametrization 
while the Gogny DIP predicts a relative change < 1%. Note that the main contribution to the 
mean free paths comes from the fact that, as shown in Fig.[IJ the Skyrme parametrization predicts 
a larger change in the absolute and relative values of the two effective masses of the spin polarized 
components. The Landau parameters and the energetic contribution of the magnetic perturbation 
[17] ] show a minor contribution to the structure functions, that in turn determine the mean free 
paths. It is worth mentioning that the Lindhard function, x^ S \ nas a rich structure in (00, q) that 
has been studied in [18j |. however, the smallnes of the magnetic perturbation is washed out in 
the response of the system by the influence of the magnetization [lfl and density effects in the 
neutron effective mass. As we can see from Fig. (3J this result shows not only quantitative but also 
qualitative differences in the neutrino transparency of magnetized neutron matter. 
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FIG. 1: Effective neutron mass at saturation density 
Po as a function of the logarithm of the magnetic held 
strength for the Skyrme SLy7 (a) and Gogny DIP 
(b) parametrizations. For each model upper (lower) 
curves correspond to spin up (down) polarized parti- 
cles. 

IV. CONCLUSIONS 

We have investigated for the first time in the context of the Landau Theory of normal Fermi 
Liquids, the effect of a strong magnetic field on the variation of the neutrino mean free path in a 
partially magnetized pure neutron system within the framework of the non-relativistic Hartree-Fock 
approximation comparing Skyrme and Gogny forces. We find that for fields up to the maximum 
strength studied in this work, B = 10 18 G, Skyrme forces show at high density an enhancement of 
neutrino transparency of the system, while Gogny forces predict a small decrease. These results 
can be explanined due to the fact that for the density, and B field range considered in this work 
the variation of Landau parameters is a minor contribution compared to the effective mass and 
magnetization. 
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FIG. 2: Relative change ratio of neutrino mean free paths for B = 5 x 10 17 G with respect to the field-free 
case as a function of density calculated with Skyrme (SLy4 and SLy7) and Gogny (DIP and D1S) forces for 
a neutrino energy E v = 15 MeV. 
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FIG. 3: Ratio of change of neutrino mean free paths as a function of the logarithm of magnetic field strength 
with Skyrme (SLy4 and SLy7) and Gogny (DIP and D1S) parametrizations at p = 3po and a neutrino energy 
E v = 15 MeV. 
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